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" l ' In a previous paper (J. Zhang et ai, J. Chem. Phys. 110, 5318 (1999)) we introduced a model 

, for polydisperse hard sphere mixtures that is able to adjust its particle-size distribution. Here we 

' give the explanation of the questions that arose in the previous description and present a consistent 

^SJ I theory of the phase transition in this system, based on the Percus-Yevick equation of state. The 

> I ■ transition is continuous, and like Bose-Einstein condensation a macroscopic aggregate is formed due 

' to the microscopic interactions. A BMCSL-like treatment leads to the same conclusion with slightly 

^5 , more accurate predictions. 
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I. INTRODUCTION 

Monodisperse colloidal suspensions are rare in nature. 
Whether these colloidal particles are artificially prepared 
or not, in the best case the particle-size distribution is 
a narrow distribution around an average size. It is ob- 
vious that such a polydisperse nature will influence the 
physics and properties of these systems, and one could 
imagine that by the use of truly asymmetric mixtures it 
is possible to create systems with a behavior that cannot 
be described by the simple monodisperse- like approxima- 
tion. 

Conceptually we can distinguish between two types of 
polydispersity. One of them, which we could refer to 
as 'intrinsic polydispersity', arises from the fact that the 
particles present in the system are different by construc- 
tion (in size, charge, or any other feature) and their char- 
acteristics are not changed by the interaction with other 
particles. This kind of systems are like multicomponent 
mixtures in which, at least in principle, the composition 
can be externally imposed. The new phenomenology 
that we can expect from this systems has its origiruip. 
fractionations into phases with different compositionsETu 
(constrained by particle conservation) and their coupling 
with other tcajisitions already present in the monodis- 
perse systemDEl. 

The second kind of jpolydispersity can be found in 
self-assembling systemaj (surfactants forming micelles, 
monomers forming chains, vesicles, etc). The aggregates 
present in these systems can be identified as the parti- 
cles, each with different size, shape, conformation, etc; 
the difference with the intrinsically polydisperse systems 
being that the composition is determined by the chemical 
equilibrium between the constituents of the aggregates. 
As a consequence the equilibrium is not constrained by 
conservation of the number of particles and therefore no 
fractionation is to be expected. In principle the system 
can compensate losses of entropy by adjusting its com- 
position. There are, however, other constraints in the 
system (the number of small constituents, for instance) 
and these may induce new kinds of transitions charac- 
terized by the appearance of one or a few macroscopic 
aggregates. As such can be cataloged phenomena as the 



^ of lamellar and columnar phases in surfac^^ 



appear aij 

solutiongj'u, emulsification failure in micro-eni«lsionsE£lc 
or long chain formation in polymer solutionstj. 

A study of polydisperse mixtures of either type is far 
from trivial. In the case of intrinsic polydispersity there 
is the experimental problem of how to fabricate colloidal 
particles according to a given particle-size distribution. 
Although in simulations this seems to be somewhat bet- 
ter under control, one easily could run into the problem 
of finite size effects due to an insufficient or inadequate 
sampling of particle sizes. This is not the case in self- 
assembling systems, although experimentally, their poly- 
dispersity cannot be easily characterized. 

Theoretical descriptions of polydisperse systems are 
mainly based on_ajSmaIl set of moments of the particle- 
size distributioroQlljiij and it is assumed that mixtures 
with the same set of moments show a corresponding be- 
havior. This seems to be a rather successful approach, 
although it is not obvious that could still be applied to 
very asymmetric mixtures. 

In two previous articles we analyzed the behavior of 
idealized versions of polydisperse mixtures. The sys- 
tem was assumed to be composed of N spherical aggre- 
gates which only interact via hard-core repulsion. The 
chemical equilibrium of the underlying constituents is ac- 
counted for by allowing particles to exchange size in such 
a way A&t the total volumell3 (compact aggregates) or 
surfacetJ (surfactant micellar membranes) of the parti- 
cles remains fixed at all times and the number of particles 
is constant. 

Under the influence of the applied pressure the 
particle-size distribution of these systems changes. In the 
case of the constrained surface, the particle-size distribu- 
tion changed from an exponentially decaying function in 
the low density limit to a single peaked distribution in 
the denser liquid phase. For sufficiently high pressures 
this system can form a number of different mechanically 
stable crystals, a monodisperse face-centered-cubic crys- 
tal, or a bidisperse AB or AB2 crystal. 

For the system with the constraint on the total volume 
of the particles, it was found by theory and simulations 
that at a rather low volume fraction of approximately 
0.26 the system could no longer be described by Percus- 
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Yevick type of equation of state, due to the formation of 
macroscopically large particles. In the present work we 
will show that this phenomenon is a true phase transi- 
tion and provide a self-consistent theory for this sort of 
system. The nature of this transition has recently been 
studied in ideal systemsEJ and found to be connected 
to Bose-Einstein condensation, and from this work we 
can conclude that interaction only changes the details, 
not the essential features. The connection of the present 
model with Bose-Einstein condensation was earlier sug- 
gested by D. FrenkelE3. 

The remainder of this paper is organized as follows. In 
Sec. H we will show the results of computer simulations 
we have performed, and explain some of the details in- 
volved. These results are compared with a theory based 
on the Percus- Yevick equation of state in Sec. III. As 
we show in Sec. |^ the treatment based on the heuris- 
tic BMCSL equation of state, leads to a slightly better 
theoretical description. In Sec. ^ we finish with a brief 
discussion on some of the issues raised in this paper. 



II. SIMULATIONS 

The system under consideration here is one formed by 
N spherical particles with different sizes. These parti- 
cles are only interacting via a hard-core repulsion, and 
hence temperature is not a relevant parameter. We allow, 
however, pairs of particles to exchange volume under the 
constraint that the total volume occupied by all particles 
remains fixed. As a consequence the system has the free- 
dom to explore different size-distribution functions-And 
relax to the "optimal" one. In the previous workli3 we 
found by computer simulations that beyond a relative 
small pressure or density, some of the particles became 
macroscopically large. From theoretical arguments an 
upper bound to the pressure was obtained for fixed vol- 
ume fraction, above which a single size distribution is 
unstable. However, simulations indicated a discrepancy 
with this result. 

The formation of macroscopic particles (in simulations 
of th|e-|Size of the simulation volume) beyond a critical 
valueuJ is in fact comparable to Bose-Einstein condensa- 
tion, or alternatively as a sort of sedimentation. In order 
to describe and simulate this phenomenon in our system 
in a correct manner, one has to eliminate the finite size 
effects that originate from the relative small number (typ- 
ically of the order 1000) of particles that is being used 
in computer simulations. This could easily be solved by 
increasing the number of particles by one or several or- 
ders of magnitude, but this is seldom a preferred option. 
Moreover, in this particular case there exists a much more 
elegant solution to the problem. 

The introduction of huge particles in a sea of smaller 
ones has two effects. The first is the reduction of the 
available volume that can be occupied by the smaller 
ones, and the second is the introduction of a hard wall 



formed by the big particles. It turns out that it is the 
second effect that caused the discrepancy between the 
simulation results and the upper bound of the equation 
of stateO. The natural solution is therefore to eliminate 
this particular effect, which can be done by extracting 
the big particle(s) out of the simulation box. 

We will assume that there is only a single big particle 
(aggregate) with a volume Vq > 0, but it can easily be ex- 
tended to include several aggregates, albeit their number 
with respect to the simulation is not relevant. The total 
volume that can be occupied by the N particles plus the 
aggregate, is denoted by Vf = Vo + N^{a^) = N^aj, 
where {a^) is the average diameter cubed of the particles 
and (7/ defines a natural length scale. The choice of ct/ 
in this formulation is based on the infinite dilute system 
for which the size of the aggregate is zero. The volume 
that is accessible to the particles is denoted by V. But 
since we have extracted the aggregate out of the system 
the total thermodynamical volume of the system is given 
by Vt ^V + Va. 

On this system we have performed Monte Carlo 
(MC) simulations using the isothermal-isobaric constant 
(NPT) ensemble. For a detailed descriptioii-on Monte 
Carlo simulations we refer the reader to Ref.EJ; here we 
will only list the main features. As positional order will 
play no role, we will assume the simulation volume V 
to be cubic and we will be using four different types of 
moves. 

The first and simplest move is randomly selecting a 
particle and displacing it isotroptically. The second type 
of move was introduced beforetZI. It selects randomly a 
pair of particles and attempts to exchange a finite amount 
of volume, by adding this to one, and subtracting it 
from the other particle. Hereafter their new volumes are 
used in order to obtain the appropriate diameters. Obvi- 
ously, if this process led to a negative volume, the move 
could not be accepted. Apart from that, both moves 
are accepted provided that no overlap is caused by these 
changes. Moreover, the amount of displacements and vol- 
ume exchanges are drawn homogeneously from an inter- 
val, of which the size is continously adjusted in order to 
obtain an average acceptance per move between 35-50%. 

In the third move we allow the simulation box to shrink 
or expand isotropically, in order for the system to equi- 
librate with respect to the applied pressure. This is in 
principle also a standard move in MC simulations using 
the NPT ensemble. However, in this case it is required 
to modify the acceptance criterion of this move. This 
is caused by the fact that the thermodynamical volume 
Vt and the volume of the simulation box y, although 
related, are not identical. 

The weight of any non-overlapping configuration in the 
iVPT-ensemble is proportional to exp(— /3PVt), where 
(3 = \/{kBT) the inverse temperature. In addition to 
that, there is in MC-simulations a correction factor 
due to the use of scaled coordinates in order to facilitate 
the volume moves and periodic boundary conditions. As 
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a result the volume move, if no overlaps are produced, is 
accepted with the probability 

Min(l,exp[/3P(t4"^ - V^°'>) - iVlog(y("V^^°^)]), (1) 

where the labels n and o refer to the new and old config- 
uration respectively. 

The fourth type of move is the exchange of volume be- 
tween a particle and the aggregate. Although this seems 
to be equivalent to the volume exchange between parti- 
cles, there is a subtle difference. Provided that no overlap 
is produced and both volumes remain positive this move 
is allowed. However, if the aggregate changes its volume, 
this results in a change of the available volume V for 
the A'^ particles, while the thermodynamic volume Vt re- 
mains fixed. As a result this move needs to be accepted 
with probability 

Min(l,exp[-Aflog(y("V"^'^''')])- (2) 

We have performed a series of simulations on this sys- 
tem with A^ — 1000 particles. In each simulation run we 
used 10^ sweeps, where in each sweep on average we try A^ 
particle moves, N/2 volume exchanges, 1 volume move, 
and 5 exchanges of volume between the aggregate and 
one of the normal particles. Average quantities where 
obtained from a single run after equilibrium was reached. 
The time required to reach equilibrium strongly depends 
on the initial configuration. In simulations where an equi- 
librium configuration of a slightly higher or lower applied 
pressure is used, usually one run is sufficient. If started 
from a monodisperse simple cubic lattice, one or several 
runs are required, depending on the initial choice of the 
size of the aggregate and volume of the simulation box. 



The results of the equation of state are shown in Fig. 
where the the reduced pressure P* ~ PPa^j is plotted as 
a function of the volume fraction 77 = V/ /Vt- The results 
obtained from expansion and compression runs, as well 
as started from monodisperse systems, all lead, within 
the estimated errors, to the same results. Near the vpl- 
ume fraction 0.26, where the critical value was foundO, 
a change in slope can be observed. 

A natural order parameter for this system is the rela- 
tive amount of fixed total volume of the particles found 
in the aggregate V^/Vf and is shown in Fig. ^j. Compar- 
ison with simulations on systems with less particles and 
the absence of hysteresis in the equation of state suggests 
the possibility of a continuous phase transition. 
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FIG. 2. The solid line is the theoretical prediction 
of the order parameter Vo/Vf. The circles are results of 
MC-simuIations on a 1000 particles system. The inset shows 
an enlargement of the area near the phase transition and in 
addition simulation results of a 512 particle system are indi- 
cated by squares. 
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■ ■■ Eq.(3) 
- Eq. (16) 
o N=1000 




FIG. 1. Equation of state. The solid and dotted line are 
the theoretical predictions above and below the transition re- 
spectively. The points are obtained from MC simulations on 
a 1000 particle system. The inset shows en enlargement of 
the transition point. 



The simulation results below the critical value are con- 
sistent with those obtained previouslylHl. Above this 
value the results are different. The presence of the aggre- 
gate seems to prevent the formation of other macroscopic 
particles. Note that the difference is not so much the fact 
that there is a aggregate, but the way it is treated as an 
external buffer, thus eliminating the surface effects. 



III. FERGUS- YEVICK THEORY 

In order to describe the system under consideration 
theoretically, one should realize that the presence of a 
thermodynamic volume Vt as well as an available volume 
V for the particles should carefully be taken into account. 
The system can be described by a polydisperse system of 
N hard spheres in a volume V , which is in contact with 
an external bath, because in addition to the exchange of 
volume between the particles themselves, this can also be 
done with this bath or aggregate. 
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We will describe the polydisperse system within the 
Percus-Yevick papproximation of a polydisperse hard- 
sphere mixturecS, from which the equation of state yields 



e'^ 1-6 (1-6)2 (1-6)3' 



(3) 



where 6 is the fcth moment of the particle-size distribu- 
tion function 



(4) 



Here p = N /V is the number density defined with re- 
spect to the available volume, rj = Vf/VT the volume 
fraction, and W[v) is the continuous particle-size distri- 
bution function, which due to the type of interaction, is 
a function of the particle volume v. The Helmholtz free 
energy functional F{N, Vt,T) of this system, however, is 
a function of the thermodynamic volume Vt 



N 



log(pA3) - 1 + / dvW{v) log W{v) + 



(5) 



where the excess free energy per particle is given by 

,2- (6) 



= -iog(i -6) + 
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6(1-6) 26(1-6)2 



In order to obtain the equilibrium distribution, W{v), 
we need to minimize this free energy functional, taking 
into account that there are two constraints which need to 
be satisfied. This is solved by adding the following two 
terms to the free energy functional 

- Cf, j dvW(v)^ Czl^j dvW{v)v + ^^ , (7) 

where £o ^-nd £3 are Lagrange multipliers. The first 
term is due to the requirement that the particle-size dis- 
tribution function should be normalized to unity, and the 
second is that the combined volume of all particles and 
the aggregate is fixed to be Vf. 

The equilibrium particle-size distribution can now be 
obtained by minimizing the free energy functional with 
respect to W{v) and the free parameter Vq. Stationarity 
with respect to changes in W{v) leads to 



log(W^(«))-log(l-6) + 
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1 - 6 



(8) 
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0-2 + -l3PcT^ -Co- Cs-a^ = 0, 
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from which we can immediately derive the functional 
form of the function W(v) that will minimize the free 
energy 



W{v) — exp 



(9) 



where the coefficients ai are determined by 
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(11) 



The values of 013 and ao are fixed by the constraint on the 
total combined volume of the particles and the aggregate, 
and the normalization of the particle size-distribution, re- 
spectively. Formally, however, the coefficient is given 

by 



a3 



-g(/3P-^3). 



(12) 



In addition to these equations we have to minimize 
with respect to the volume of the aggregate. Note that 
the 6 depend on the value of Vo through the number 
density. Therefore we obtain 



dl3F d(3F 



dVo 



dVr 



-£3^I3P-C3 



0. 



(13) 



This result does not automatically imply that a3 = 0. 
The reason is that we have to minimize the free energy 
with respect to the volume of the aggregate under the ad- 
ditional condition that the aggregate has a non-negative 
volume. As a consequence this can result in the possibil- 
ity that the minimum is found for Vq = and the sta- 
tionarity equation (|l^) not being satisfied. This leads to 
a natural splitting of the behavior of this system. For low 
pressures the volume of the aggregate will be zero, and 
the particle size-distribution is the exponential of a third 
order polynomial in the particle diameter. On increasing 
the pressure the coefficient as, which can only be non- 
positive in order to allow a proper normalization, goeSpto 
zero. In this region the behavior is as explained befordlZI. 
According to the theoretical description in the Percus- 
Yevick approximation, the coefficient as becomes zero at 
volume fraction t] = 0.260198 pressure P* = 1.343442. 
At this point the value Vo = not only minimizes the 
free energy, but also satisfies the stationarity equation 
( p^ ) . For pressures larger than this value the stationarity 
equation can always be satisfied with a positive volume 
for the aggregate and hence the coefficient 013 remains 
zero. 

It turns out that this point is the location of a con- 
tinuous transition. Beyond this point the system relaxes 
by transferring volume to the aggregate, and effectively 
scaling the polydisperse system to a smaller size. If we 
denote by a script c the values of the critical system, this 
can be shown in the following way. If we replace the 
coefficients ak in the particle size-distribution by q'^ack, 
where g is a positive number, and fix ao by the normal- 
ization to unity, it follows directly from the special form 
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of W{v) that the moments of the size-distribution are 
given by {a'') = (<j^)cl(t ■ This resuhs in 



ax 



„3 3r;, (g^) 

and provided that q is determined by 

3 r/ 1 - 
9 = :j 



(14) 



(15) 



this is a sohition of Eq. (|T^. Applying this procedure 
to ai leads to the same expression for q. Therefore the 
scaled size-distribution function leads to a minimum free 
energy, albeit for a different applied pressure. 

The pressure equation of state beyond the critical point 
also becomes fairly simple now. If we rewrite the pressure 
(^ using the solutions ( p^ and ([Tl|), we obtain 



2r] 



(16) 



where we used that if 0:3 = the functional form of W{v) 
leads to the identity 3{a°) + ai{a) + 2a2((7'^) = 0. We can 
also derive a simple expression for the order parameter 
Vo/Vf of the system 



Vf 



= 1 - ^. 



(17) 



The curves of the equation of state and the order param- 
eter are shown in Figs. |l| and ^ respectively, and show a 
perfect agreement with the simulation results. In Fig. |l| 
the dotted line represents the equation of state below 
the transition, i.e. Eq. (H) , while the solid line is the one 
above as given by Eq. (^6|) . The later one is extended to 
the zero density limit in order to visualize the difference 
between both branches. 

From these results two additional illustrative proper- 
ties can be derived for the system beyond the point of 
the transition. Using that the particle size-distribution 
only changes according to a dimensional scaling or from 
the order parameter (O), we find that 



(18) 



This simply states that the system becomes scale invari- 
ant, with a constant pressure if it is defined in terms of 
the average particle volume. The other result is that ^3, 
which can be identified with the local volume fraction of 
the particles, remains fixed, i.e. ^3 = rjc. These results 
can be used in order to determine the point of transi- 
tion from simulations. This is illustrated in Fig. ^, where 
P*{a^)/aj as a function of rj, and the applied pressure 
P* as a function of ^3 are shown, as measured in our MC- 
simulations. According to these simulations the transi- 
tion is found at a volume fraction rj = 0.2623 ± 0.0003 
and pressure P* = 1.360 ± 0.003. 




FIG. 3. Comparison of the simulation results for a 1000 
particles system with the theoretical predictions. The circles 
represent the applied pressure P* as function of the volume 
fraction rj. The squares the local pressure P* (a^) / as func- 
tion of r], and the diamonds the applied pressure P* as func- 
tion of the local volume fraction N{v)/V . 



The prediction of the particle size-distribution W(v) is 
compared with the simulation results in Fig.^. By plot- 
ting ln{W(v)/ P*) as a function of P*v, the distributions 
above the transition coincide and agree perfectly with the 
predicted curve from Percus-Yevick. The discrepancy at 
higher values for P*v is merely a consequence of an in- 
sufficient sampling of the largest particles in the system. 
For P* — 1.50, there is the additional effect that this is 
close to the transition, where the system tends to probe 
larger particles. 
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FIG. 4. The particles size-distributions are shown by plot- 
ting ln{W{v)/P*) as a function of P*v. Due to this scaling 
all distributions above the transition should coincide, without 
using the actual values of the pressure and volume fraction at 
the transition. The solid curve is the predicted curve by the 
Percus-Yevick approach. 
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It turns out that within the Percus-Yevick approxi- 
mation the system shows a continuous phase transition, 
which is characterized by the non-zero volume of a aggre- 
gate. In fact it shows some similarities to Bose-Einstein 
condensation, because what actually happens is that the 
interaction of microscopic piarticles leads to the formation 
of a macroscopic aggregatalS. 



The derivative of the free energy with respect to the vol- 
ume of the aggregate leads again to 



dl3F d(3F 



dVo 



dVr 



£3 = WP - 



(24) 



IV. BMCSL-THEORY 

It is a well known fact that the Percus-Yevick equa- 
tion of state in a monodisperse hard sphere system can 
be improved by combining it with the virial equation of 
state—This results in the Carnahan-Starling equation of 
statetj. A similar approach can be made for polydisperse 
sphere mixtures, and the resulting heuristic equation of 
state is due to BoubL fk . M ansoori, Carnahan, Starling, 
and Leland (BMCSL)EJe 
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(1-6)^' 

(19) 



A straightforward integration of this equation of state 
shows thajt the excess free energy up to a constant is 
given byEHl 



64' 



366 



1 ) ln(l - 6) 

6^ 



' 6(1-6) ' 66(1-6)^"' ^^^^ 

The additional constant is independent of the volume, 
but could in principle depend on the moments of the 
size-distribution, without changing the equation of state. 
Such a constant would be relevant for this particular 
model, since the size-distribution of the system can be 
changed, and therefore an additional change in free en- 
ergy can be produced. From general applicability of the 
free energy functional to different models, however, it can 
be shown that this contribution necessarily vanishes. 

We now can apply the same procedure as was used in 
the case of the Percus-Yevick approximation. The func- 
tional derivative of the free energy with respect to W{v) 
leads to the same functional form (^ as before, but the 
stationarity equations on the coefficients ai need to be 
modified 

36 



ai 



0L2 



(1-6 
3^2 

0^^6(1-6)2 
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(21) 



log(l-6), (22) 

^3 



3g 
■2 
^3 
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(/3P - l^z) 



21og(l-6 



6(2-6 



(1-6) 



(23) 



The second term on the right hand side of Eq. (gSj) is, 
except in the zero density limit, always positive. As a 
consequence Eq. can never be satisfied, because such 
a solution would lead to a positive 0^3 and hence a non- 
normalizable size-distribution function. The minimum 
free energy for given volume fraction will therefore not 
be a solution of all the stationarity equations simultane- 
ously. 
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FIG. 5. The Helmholtz free energy per volume as a func- 
tion of the volume fraction for fixed order parameter Vo (Vf , 
i.e. constant size of the aggregate (solid curves) and for scaled 
distribution functions (dotted curves). Both sets are similar 
but use only a different parameter for characterization. 



Since the stationarity equations (|2l| ) - (23) are directly 
related to the equations for chemical equilibrium between 
the different particle sizes, they need to be fulfilled. If 
we fixed the volume of the aggregate Vb , these equations 
would be sufficient to obtain the unique distribution func- 
tion that minimizes the free energy. In Fig. ^ we have 
plotted the free energy per volume for several fixed val- 
ues of Vb (solid curves). On these curves the stationarity 
equations are satisfied, hence the value of Vb that leads 
to the minimum free energy corresponds to the equilib- 
rium situation. These curves are bounded by an lower 
and upper limit of the volume fraction. The lower limit 
corresponds to the case where a\ — a2 — ^ and the up- 
per limit corresponds to the case where a-i = 0. Outside 
this interval no self-consistent solution of the stationarity 
equations can be obtained. 
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Similar to what we found in the Percus-Yevick ap- 
proach, we can scale the particle size-distribution func- 
tion such that we obtain a solution to the stationarity 
equations for a different volume fraction by replacing ak 
by q'^ak and fixing the normalization constant. The free 
energy per volume for several of these scaled solutions 
is also shown in Fig. || as dotted curves. Both sets of 
curves are similar, but use a different parameter to de- 
scribe them. 

The equilibrium situation of our model is the one that 
corresponds to the lowest free energy for given volume 
fraction and will correspond to the envelope of the set 
of curves in Fig. ||. One can show that we have two 
regions, which is illustrated by this figure. Below the 
transition the equilibrium situation corresponds to the 
curve for which Vq — 0, i.e. without a macroscopic ag- 
gregate. Above the transition we will follow the branch 
of a scaled solution, emerging from the end point of this 
curve determined by = 0. The transition point is 
therefore obtained by self-consistently solving the case for 
Vo = 0:3 = and leads to a volume fraction 77 = 0.262611 
and pressure P* = 1.356275. 

Actually the same reasoning could have been used for 
the Percus-Yevick approach. But since for that case 
all stationarity equations could be solved simultaneously, 
there was no need to follow this line of argumentation. 

The numerical predictions of the transition point for 
the BMCSL approach are slightly better in agreement 
with the simulation results than the results from Percus- 
Yevick. The difference between both predictions is how- 
ever too small in order for it to be visible in any of the 
figures. In principle also the size-distribution function 
can be used to distinguish both theories, however, in or- 
der to do so the number of particles should be increased 
several orders of magnitudes in order to obtain sufficient 
statistics for the larger particles and to make the differ- 
ence in the tail of the distribution visible. 



V. DISCUSSION 

We have demonstrated how to obtain a self-consistent 
theory for a polydisperse hard sphere system, where par- 
ticles are able to change their size by exchanging vol- 
ume. Within the Percus-Yevick approximation this sys- 
tem shows a continuous phase transition, in which the 
microscopic interactions lead to the formation of a macro- 
scopic aggregate, a behavior similar to Bose-Einstein con- 
densation. The theoretical results have been compared 
with a new series of computer simulations and show an 
excellent agreement, because in the new treatment we 
have eliminated the finite size effect that previously led 
to discrepancies between theory and simulation. 

This peculiar system, which at best can be considered 
to be an extremely idealized version of spherical micelles, 
exposes an unusual behavior in the BMCSL-approach. 
The equilibrium state of the system is one which does 



not satisfies the stationarity equations but is determined 
by the boundary conditions on the problem, i.e. a nor- 
malizable particle size-distribution. Apart from that the 
difference between both approximations is minor and in 
order to determine which one is better, simulations would 
be required with several orders of magnitude more par- 
ticles. 
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